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ABSTRACT 


A new method is presented for estimating the reflection of a random, multi¬ 
directional sea from a coastal structure. The technique is applicable to an array 
of wave gauges of arbitrary geometry deployed seaward of the reflector. An 
expansion for small oblique wave incidence angles is used to derive an 
approximate relationship between measured array cross-spectra and a small 
number of parameters that describe the incident wave properties and the 
reflectivity of the structure. Model tests with simulated array data demonstrate 
that for wave incidence angles less than about 30° the new technique provides 
accurate and robust estimates of the gross properties of incident and reflected 
waves. 

The new method is applied to array data acquired offshore of a permeable, 
rubble mound breakwater in Monterey Bay, California. The estimated reflection 
coefficients decrease approximately linearly with increasing frequency. Whereas 
the observed reflections depend only weakly on the incident wave energy, the 
fraction of the incident wave energy flux transmitted through the breakwater 
decreases with increasing wave energy, suggesting that dissipation is enhanced 
with large amplitude waves. 
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T1i,Tlr=surface elevation functions of incident and reflected 
waves 

k,lc,=vector waveniiinbers of incident and reflected waves 
dz=incident wave amplitude 
x=horizontal position vector 
t=time 
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<|>=phase lag between incident and reflected waves 

h^=cross-spectrum of sensors m and n 
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Ro/ <l>o/ R 2 / 4 > 2 =^^eflection parameters 
L=distance from coordinate origin to the breakwater 
E=incident wave energy spectral density 
0^„=mean propagation direction of incident waves 
0r„=root -mean-square directional spread of incident waves 






^0* $o=estiinates of incident and reflected wave 
parameters 

fiM,=observed and model cross-spectra 
e=root-mean-square misfit norm 
fi=average of the auto-spectra 
f=frequency (Hz) 

Fi=incident energy flux density 
Fr=reflected energy flux density 
Ft=transmitted energy flux density 
Fd=dissipated energy flux density 

Cgi, Cg2=group velocity outside and inside the harbor 
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I. ZNTRODDCTZON 


The interaction of ocean surface waves with a permeable, 
rubble mound breakwater is complex. In general, the 
incident waves are partially transmitted through or over the 
porous structure, partially dissipated by breaking on the 
rough slopes and by turbulent friction within the 
breakwater, auid partially reflected. The reflection of 
waves from breakwaters (and other reflectors such as sand 
bars and seawalls) is usually measured with an array of 
pressure sensors or surface height gauges deployed seaward 
of the reflector (e.g., Thornton and Calhoun 1972; Mansard 
and Funke 1980; Yokoki et al. 1992). The interpretation of 
these measurements is complicated owing to the fact that the 
incident and reflected waves are phase-coupled. 

In many studies the array analysis is simplified by 
assuming that the incident waves are \ini-directional, 
propagating perpendicular to the reflector. Thornton and 
Calhoun (1972), Morden et al. (1976) and Goda and Suzuki 
(1976) used two sensors positioned on a line perpendicular 
to the reflector to decompose the wave field into incident 
and reflected wave contributions. This technique breaks 
down at the frequency where the wavelength is equal to twice 
the sensor spacing. Meuisard and Funke (1980) overcame this 
problem by applying a least-squares-fit technique to three 
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sensors, cind this approach was extended to linear arrays 
with an arbitrary number of sensors by Zelt and Skjelbreia 
(1992) . 

The ass\imption of normally incident waves used to 
estimate reflection is often violated in a natural coastal 
environment. Refraction of swell over complex bathymetry 
may result in significemtly oblique angles of incidence at 
the breakwater, especially if the breakwater is not aligned 
with the surrounding depth contours. Directionally broad, 
locally generated seas are typically only weakly refracted 
in depths greater than 10 m and can approach a breakwater at 
relatively large oblique angles. Although in principal 
array measurements seaward of a reflector can be used to 
infer the reflection of a directionally spread wave field, a 
very large number of sensors is required to obtain reliable 
estimates of the directional spectra of both incident and 
reflected waves (Isobe and Hondo 1984). 

In this study, a new method is presented for estimating 
wave reflections from array data. The assumption of normal 
incidence used in previous studies is relaxed, but angles of 
incidence are assumed to be small. Reflection from the 
breakwater is assximed to be a linear process governed by 
Snell's law (i.e., specular reflection). An expansion for 
small oblique incidence angles is used to derive approximate 
relationships between the array cross-spectra and a small 
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number of parameters that describe the incident wave 
properties and the reflectivity of the breakwater. 

The new estimating technique is applied to pressure 
array data acquired offshore of the Monterey Harbor 
breakwater at Monterey, California. The field site, array 
geometry and data acquisition are described in Section 2. 
Section 3 presents the new estimation technique, and 
demonstrates the accuracy and robustness of the method 
through model simulations. Results of the analysis of 
reflections observed at the Monterey breakwater are given in 
Section 4, followed by a discussion and conclusions in 
Section 5. 
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ZZ. SXPBRZNBNT 


The experiment was conducted at the Monterey Harbor 
breakwater in Monterey, California (Fig. 1). Situated at 
the southern end of Monterey Bay (Fig. la), the permeable 
breakwater is exposed only to swell with deep water 
propagation directions ranging from west to north-northwcst 
cuid seas generated by local northerly winds. The 
(refracted) directional wave spectra at this shelte-_d site 
are typically narrow with most of the energy acoroaching the 
breakwater from northerly directions within abr it 20° of 
normal incidence to the breakwater. 

Construction of the 400 m long rubble mound breakwater 
was completed in 1934, and that same year a 120 m long 
extension was added (Fig. lb). During subsequent years, 
much of the interior of the breakwater washed away and the 
elevation of the structure decreased so that it was often 
overtopped at high tide (Thornton and Calhoun 1972). In 
1993, the breakwater extension was improved by increasing 
the elevation back to the original design height of +3 m 
relative to MLLW, and by adding on the harbor side a single 
cover layer of 8-ton-average amor stone. A cross section 
of the breakwater can be found in Thornton and Calhoun 
(1972). The present field study was conducted at the 
improved extension. 


4 







A 2-dimensional array of six pressure sensors was 
deployed seaward of the breakwater in a depth of about 16 m 
(Fig. lb) to measure both the incident and reflected wave 
fields. The small depth variations (less than 1.5 m) across 
the array are neglected. The array aperture is 31 m x 15.6 
m and the distance between the breakwater and the closest 
sensor P6 is about 52 m. Additionally, a single pressure 
sensor P7 was positioned inside the harbor at a distance of 
about 47 m from the breakwater to measure the transmitted 
wave field. Sensor locations relatively close to the 
breakwater were chosen to avoid contamination by edge 
effects (e.g., diffraction around the tip of the 
breakwater). 

The array data were acquired with a tattletail micro¬ 
processor located in the center of the array, which was 
ceLbled to a radiotelemetry system on the breakwater, from 
where the data were transmitted to a computer at the Naval 
Postgraduate School. Data were acquired with a 2 Hz sample 
rate continuously for eight months between August 1993 and 
March 1994. 

Typical power spectra measured along a line 
perpendicular to the breakwater, converted to surface height 
with a linear theory depth correction, are shown in Fig. 2. 
The swell spectrum is most energetic in the frequency range 
of about 0.06-0.1 Hz. Large differences in energy levels 
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between the different sensors and pronounced hills and 
valleys in the spectra indicate nodes and antinodes as 
discussed in Thornton and Calhoun (1972). At frequencies 
below about 0.05 Hz relatively energetic infragravity 
motions are apparent. Energy levels inside the harbor 
(sensor P7) are significantly reduced at swell frequencies 
owing to reflections and/or dissipation. 
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ZZZ. BSTZM&TZON OF RBFLICTZONS 


In this section, a new method is presented for 
estimating reflections of ocean surface waves from coastal 
structures (e.g., brea)cwaters, seawalls) or natural shores 
(e.g., sand bars, steep foreshores) using an array of 
pressure sensors or surface height gages located seaward of 
the reflector. The reflector profile and bathymetry are 
assijmed to be approximately uniform along the reflector and 
depth variations in the array vicinity are neglected. 

The surface elevation function Tii(x,t) of a random, 
stationary incident wave field can be expressed as a 
Fourier-Stieltjes integral: 

T],(x,t)= f f exp[i(Jt'X~<i>t) ] ciz(<i),k) (1) 

Jut J k 

where dz(0),k) is the complex valued anplitude function, (O is 
the radieui frequency, k is the vector wavenumber, x is the 
horizontal position vector (the coordinate frame is defined 
in Fig. lb) euid t is time. Assuming that reflection is a 
linear process, the surface elevation function Tir(3c,t) of 
the reflected waves can be expressed as: 

Tj,(x,t)=f f J? exp[i(Jfc,-x-cdt+<j>) ] d2{a),l:) (2) 

J» JJk 

where R is the reflection coefficient, <|> is the phase lag 
between incident and reflected waves relative to x= 0 , euid kr 
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is the vector waveniimber of the reflected wave. R, <]) and kj. 
are generally functions of k and to. Adding Eqs. 1 and 2 
yields the total surface elevation function Ti(x,t): 


n (x, t) = j* (*‘X-co t) ] 

+R exp [i (Jt/ x-(o t+<j)) ] ] dz{ui,k) 

The cross-spectrum h^((D) of two sensors at positions 
and is defined as: 

f exp[i«T]h^(o))dto s £{ii (Ji^, t) t}* (x., t+t)} 

where E{} denotes the expected value and * the complex 
conjugate. Substitution of Eq. 3 in Eq. 4 yields: 

“ /Je3cp [iibxj +R exp [i (k^- I - 
X [exp[-iJrxJ +E exp[-i (Jt,* Xa+<|>) ]]S{v>,k) dk 

where 8 ( 00 ,k) is the waventimber-frequency spectrum: 

Sim k) - dz*((o,k)} 

' diidk 

The wavenumber magnitudes of both incident and reflected 
waves are given by the linear dispersion relation 
CO*=gk[tanh(kh) ], where k=(k |=jk,j, g is gravity and h the 
water depth. Introducing polar coordinates 
(k=: [k(cos (6)) ,k(sin(6)) ]) and further assuming that the 
reflection is specular (kr=[-k(cos(0)),k(sin(0))1), Eq. 5 


(5) 


( 6 ) 


8 



can be expressed as (dropping the frequency dependence): 


f ^0^(6) Sid) dd (7) 

" 2 

where 0^(0) is given by: 

=lexp [iic{x;„cose+y^ine) ] 

+i?(0) exp [i (-icXj,cosO+icy„sin0+(l) (0)) ] ] 

X [exp (-iic (x;,cos0+y„sin0) ] ^ ®) 

+i2{0) exp [-i (-Jcx^cosO+icy„sinO+<|){0)) ]] 

$ 

auid S(0) is the directional spectrvun of incident waves. 

Asstuning that incident waves approach the breakwater at 
small oblique euagles, R and ^ can be expanded for small 0: 

R{Q) =Ro+ie202+22^e*+. . . (9) 

<|)(0) =<|)o+<|)202+4)^0«+. . . (10) 

The odd-order terms in these expansions vanish because R(0) 
and <j»(0) are symmetric fxinctions. Sxibstituting Eqs. 9 and 
10 in Eq. 8 euad expanding G^(0) for small 0 yields : 

with the lowest three coefficients A^, B,„ emd given by: 


A^= exp [iic(x„-x„) ] +2RoCOS [k{x„+x„) -^q] 
+i2o^exp [ -ik{x„-x„) ] 


(12a) 
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(12b) 


Cn«=—(y«,-yn) )exp [iA:(x„-x„) ] 


+ [ 2 / 22 (y*-yn) ^3 cos(ic:{x„+x„) -<|)o) 

+ [ 2 J?o<|» 2 +i 2 o^(-Jf*+-^n) 3 sin(Jc(x„+x„) -4)o) 

+ [-|j2o*(-ic^ (y„-y„)^+ik{x„-x„) )+2i?o/?2] exp [-iJc{x„-x„) ] 


(12c) 


Gna(0) caul be accurately approximated by a truncated 
expansion of the form Eg. 11 if all the expeuision 
coefficients (Eg. 12) are jSOd) . Thus the technigue 
presented here is generally valid only for contact arrays 
(i.e., |k • x|^0(l)) positioned close to the reflector 
(i.e,, within a wavelength so that <j>:SO(l)). 

The exact Go„(0) (Eg. 8) of sensor pair n=3, m=l (Fig. 

lb) for f=—=.08 Hz (the domineuit swell freguency) is 
2is 

con^ared to the linear (A„a+Bna0) and guadratic (A™+Bn„0+C^0*) 
approximations in Fig. 3. In this calculation the 
brealcwater is idealized as a partially absorbing wall at 
x=L=83 m with R independent of 0 euid (|>=2kL[cos(0) ] (Fig. 
lb). The linear approximation diverges from the exact 
solution for |0|>5-15° but errors are generally within 10% 
even for incidence angles as large as 30°. The guadratic 
approximation is more accurate for small values of 0, but 
diverges sharply from the exact solution for large 0 with 
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deviations exceeding 20% for |0|S3O°. Similar calculations 
for other sensor pairs and other frequencies in the swell 
range (0.05-0.12 Hz) show comparcdjle agreement between the 
exact Goal(6) smd the linear and quadratic approximations. 

Substitution of the trxmcated expansion of G^(6) (Eqs. 
11,12) in Eq. 7 yields the quadratic approximation of the 
cross-spectrum h„a,; 

where E, 9-. ,.., and 0^ are the energy spectral density 
(integrated over all directions), mean propagation direction 
and root-me 2 ui-square directional spread of the incident 
waves: 



0 , 


mean 



(14a) 


(14b) 
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(14c) 
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Eq. 13 relates the array cross-spectra h^, to three 
parameters of the incident wave field (E, 0rma) and 

four reflection parameters (Ro , <l>o * ^2 » <l >2 ) • general, 
inverse algorithms C 2 Ln be developed that search for a 
combination of these seven parameters on a frequency by 
frequency band basis such that the associated cross-spectra 
are as close as possible to the observed cross-spectra, but 
this requires 2 ui extensive array. Since the six-element 
array deployed in the present study is relatively small in 
aperture, the contributions of the quadratic (Cna.) terms to 
Eq. 13 is small even for wave incidence cungles as large as 
20° (e.g., coitpare the linear and quadratic approximations 
of G3 i( 0) in Fig. 3). Thus, the higher-order parameters 0^, 
Rj, and <j>j may not be resolvable within the uncertainty of 
the array cross-spectra. Neglecting the quadratic terms in 
Eq. 13 reduces the inverse problem to only four unknowns; 
two incident wave parameters (E, 0. ^^.. ) and two reflection 
parameters (Rq , <|>o) • Estimates t, fto an^ $o were 

obtained by minimizing a single root-mean-square misfit norm 
e: 




(15) 


\ H fli / 

with the observed cross-spectra and the linear model 
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predictions: 

Since both the number of sensors emd the number of unknovms 
is small, this minimum misfit can be evaluated by 
essentially sweeping through the entire parameter space. 

For all possible combinations of fi, &<, and $ 0 / the 

cross-spectra and the misfit e were calculated with Eqs. 
12a, b, 15 and 16 to obtain a global minimum for e. The 
range of physically plausible values for fi, emd $0 

used in these calculations is 0.25E-4E, -30°-30°, 0-1, and 
0-360°, respectively, with E the average of the auto spectra 
fioa* minimum value of e was calculated by sweeping 

through all possible combinations of E, 

stepping with increments of 0.15E, 5°, 0.1, emd 15°, that 
were somewhat coarse owing to limited computing resources. 
The accuracy of the solution fi, fto' $0 was improved by 

sweeping with smaller step sizes (0.03E, 2°, 0.02, 3°) 
through a reduced parameter range (fi±0.3E, ^... ,,^^ ±10°, Ro±0.2, 

$o±30°) . 

The accuracy of the new estimation technique was 
verified through a series of model tests with simulated 
array cross-spectra. For a chosen incident swell 






directional spectrum of the form: 

S{Q )« cos^°° 

the true cross-spectra h^ of the Monterey array (Fig. lb) 
were evaluated with Eqs. 7 cuid 8. As before, the breakwater 
was idealized in these tests as a partially absorbing wall 
at x=L=83 m with R independent of 6 and (|>=2kL[cos (6) ] . In 
each test, five reindom realizations of cross-spectra finm 
were generated using the procedure described in Long and 
Hasselmeuin (1979) . Errors in the simulated include both 
statistical uncertainty resulting from finite length data 
records (160 degrees of freedom) and uncorrelated instrument 
noise (noise to signal ratio 0.05) and are roughly 
representative of the actual array measurements presented in 
Section 4. The simulated figa were then treated in exactly 
the same way as cross-spectra obtained from ocean 
observations, and the minimum misfit e and the optimal 
parameters were obtained with the inverse 

algorithm described above. 

Results of model tests for typical swell with frequency 
f=0.086 Hz are presented in Figs. 4 and 5. Similar 
cos^eurisons (not shown) for other frequencies in the swell 
band (0.05-0.12 Hz) yielded conqparable agreement. Estimated 
(ft) and true (R) reflection coefficients are in good 









agreement for both normally (Fig. 4a) and obliquely (Fig. 

4b) incident swell (discrepauicies generally within ±0.05). 
Estimates of the incident wave energy E (Fig. 4c) and 
propagation direction 9 ,^., (Fig. 5) also agree well with the 
input spectrum (errors generally less than 20% and 4°, 
respectively). The values of the misfit e, a normalized rms 
measure of the discrepancies between the model (fi^j,) and 
"observed" (fi^) cross-spectra (Eq. 15), range from 0.02 to 
0.08 (Fig. 4d). The misfits are generally larger for the 
simulations with obliquely propagating waves than for 

normally incident waves (e.g., compare the e for 6 _ =0° and 

20° in Fig. 4d) owing to errors in based on a small 0 
approximation (Eq. 16). However, aside from a slight bias 
in ft (Fig. 4b), the estimates of incident auid reflected wave 
parameters do not appear to be significamtly degraded for 
0-.«=2O° (Figs. 4,5). 

Overall, the model simulations demonstrate that the 
estimation technique is relatively insensitive to errors in 
the data and the model, and can extract accurate estimates 
of wave reflections from the array data acquired in the 
present study. 






IV. OBSERVED REFLECTIONS 


The reflection estimation technique described in the 
previous section was applied to the array measurements 
acquired at the Monterey Harbor breakwater (Fig. 1). Nine 
days were selected for analysis, that approximately span the 
ramge of conditions encountered during the experiment. For 
each of these nine days, two three hour data runs were 
processed, one acquired at low tide and one at high tide, 
with the objective to examine the sea level dependence of 
breakwater reflections. Cross-spectra with a frequency 
resolution of 0.0078 Hz and 160 degrees of freedom were 
computed for each of the 18 data runs. For every frequency 
band in the dominant swell range, 0,05 Hz-0.12 Hz, the 
inverse algorithm was applied to the cross-spectra to 
estimate the incident wave spectral density E(f), the mean 
incident wave propagation direction (f), the reflection 

coefficient Ro(f) and the phase lag between incident suid 
reflected waves * Frequencies less than 0.05 Hz and 

greater thcui 0.12 Hz are not considered here because 
infragravity waves dominate the spectra below 0.05 Hz (e.g.. 
Fig. 2; Okihiro et al. 1992) and above 0.12 Hz (i.e., 
frequencies greater than about twice the spectral peak 
frequency) local nonlinear effects may be significemt. 
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The variability of incident wave conditions was small. 
Estimates of the incident swell variance ranged from 3.7 cm^ 
to 51 cm^ (i.e., significant wave heights of 7-29 cm). 

These low energy conditions are typical for this site owing 
to sheltering effects (Fig. 1), and the fact that no major 
storm occurred during the eight month data acquisition 
period. Estimates of the mean frequency and propagation 
direction of incident waves (averages of f cund (f) over 
the swell band, weighted by E(f)) ranged from 0.065-0.089 Hz 
and from 0-18°. These incidence angles are well within the 
r 2 mge for which the present technique (based on a small 6 
expansion) is expected to be accurate (Figs. 3-5). 

Estimates of Ro(f) obtained from different data runs are 
remarkably similar (Fig. 6), even though the incident wave 
spectral levels varied by more than an order of magnitude. 
Reflection of small amplitude swell from the breakwater is 
apparently insensitive to the incident wave amplitude, 
consistent with the assumption that reflection is a linear 
process. In all cases the reflection estimates show a 
strong frequency dependence with Ro(f) decreasing 
approximately linearly with increasing frequency from about 
0.7-0.8 for f=0.05 Hz to 0.2-0.3 for f=0.12 Hz. In contrast 
to estimates reported by Thornton and Calhoun (1972), the 
present observations do not suggest a strong dependence of R 
on tidal sea level variations. 
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The misfit e between and (Eq. 15) is shown in 
Fig. 7 as a function of frequency. In the most energetic 
part of the spectrum (0.06-0.1 Hz), e is approximately 0.03- 
0.07, comparable to the misfits obtained in model tests. At 
frequencies below 0.06 Hz and above 0.1 Hz, where energy 
levels are relatively low, the misfits are slightly larger 
(0.05-0.16) than expected from model tests, possibly owing 
to directional spreading and/or nonlinear effects. 

Estimates of the phase lag <|>o(f) between incident and 
reflected waves obtained from 18 different data runs are 
approximately equal, increasing with increasing frequency as 
expected from theory (Fig. 8). Neglecting depth variations 
(i.e., changes in k) seaward of the breakwater, the 
theoretical phase lag <j>o for small incidence angles is equal 
to 2kL. This crude approximation of <|>o, taking L to be the 
distance to the crest of the breakwater, is in good 
agreement with the estimates across the entire swell band 
(Fig. 8). 
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V. DISCUSSION AND CONCLUSIONS 


A single sensor (P7) was deployed inside the harbor 
(Fig. lb) to obtain crude estimates of wave transmission 
through the breakwater. Assiiming that the propagation 
directions of incident, reflected and transmitted waves are 
nearly perpendicular to the breakwater, the fluxes of 
incident (Fi(f)), reflected {Fr(f)) and transmitted (Ft(f)) 
energy are approximately given by: 

Fi{f)^E{f)Cgj^{f) (17a) 


F^{f)=R^{f)E{f)Cg^{f) (17b) 

F,(f)=Ep,(f)Cg 2 (f) (17c) 

where Cgi and Cgj are the group velocities at the offshore 
array euid sensor P7, E is the estimated incident wave 
spectrum and Ep, is the transmitted wave spectrum measured 
by sensor P7. The transmission coefficient T(f), defined 
as: 



is plotted as a ftinction of frequency in Fig. 9 for all 18 
data rtans. The observed transmission coefficients are more 
varied^le than the reflection coefficients (Fig. 6), but do 
not suggest a consistent sea level or frequency dependence. 
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Estimates of the residual energy flxox : 

F^{f)=Fi{f)-F^{f)-F,{f) (19) 

that is dissipated through wave breaking on the rough 
breakwater slope and/or turbulent friction inside the porous 
structure, are variable, but F^/Fi generally increases with 
increasing frequency (Fig. 10). 

Estimates of the bulk incident, reflected, transmitted 
and dissipated energy fluxes, obtained by integrating Fi(f), 
Fr(f), Fc(f) and Fd(f) over the entire swell band (0.05-0.12 
Hz), are presented in Fig. 11. The observed fraction of the 
incident energy fliox that is reflected from the brealcwater 
(Fig. 11a) varies between cibout .2 and .5. These changes in 
the breakwater reflectivity are primarily the result of 
variations in the dominant swell frequency (0.06-0.09 Hz, 
Fig. 6). On the other hand transmission of wave energy 
through the breakwater appears to be a strong fianction of 
the incident energy flux (Fig. 11b). On days with very low 
anplitude swell about 40-60% of the incident energy fliix is 
transmitted through the breakwater and dissipation is weak 
(0-40%, Fig. 11c). On more energetic days about 40-60% of 
the incident energy flvix is dissipated (Fig. 11c) and the 
tremsmission is reduced to edDout 20-30% of the incident 
energy flux (Fig. 11b). 

Estimates of energy transmission and dissipation also 
show a dependence on sea level. At high tide the 
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transmitted (dissipated) energy fliixes are slightly larger 
(smaller) than at low tide, possibly owing to the fact that 
the effective width of the breakwater near the sea surface 
is smaller at high tide than at low tide. 

Waves transmitted through the breakwater may undergo 
partial reflection at the other side of the harbor and the 
associated standing wave patterns may contribute significant 
errors to estimates of Ft(f) based on Eq. 17c. Accurate 
estimation of wave transmission requires an array of sensors 
on the harbor side of the breakwater, which was not 
availaible in this study. Furtheinnore, the range of 
conditions encovntered in the present experiment was rather 
limited, and the reflection/transmission coefficients 
observed under benign conditions may not be representative 
for the breakwater performance under storm or large 
amplitude swell conditions. More extensive measurements are 
needed to evaluate the performance of permeable rubble mound 
breakwaters. 

The main result of the present study is the development 
of a new technique for estimating the reflection of a 
random, directionally spread wave field from a coastal 
structure (e.g., a breakwater or seawall) or natural sand 
bars. The estimation technique can be applied to a compact 
(i.e., aperture less than a wavelength) array of pressure 
sensors or surface height gauges of arbitrary geometry. 
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deployed seaward of the reflecting surface (within a 
wavelength). Model tests demonstrate that for wave 
incidence angles less than cibout 30° the new method cem 
provide accurate estimates of the gross properties of 
incident and reflected waves. 












APPBIDZX 


Flgur* 1. a) Location of the Monterey Harbor breakwater, 
b) Coordinate frame and locations of pressure sensors 
deployed near the 120 m breakwater extension. The three 
sensors in the center (P2,P4,P5) form an equilateral 
triangle with dimension 2.5 m. PI is 15 m seaward of the 
trieuigle, P3 cuid P6 are both positioned at a distance of 13 
m from the center. Sensor P7 is inside the harbor. 
Soiindings are in meters relative to MLLW. 
rlgur* 2. Typical power spectra of array sensors PI 
(solid), P5 (dashed), P6 (dotted), positioned on a line 
perpendicular to the breakwater (Fig. lb). The dash-dot 
curve is the spectrum measured inside the harbor (sensor 
P7) . 

Figure 3. The exact Gna(0) (solid lines, Eq. 8) compared to 
quadratic (left panels) euid linear (right peuiels) 
approximations (dashed lines, Eq. 11), for sensor pair n=3, 
m=l, with the reflection coefficient, R, equal to 0, .5 and 
1 . 

Figure 4. Results of model tests for normally (0 m..» =O°) and 

obliquely (0_,=20°) incident waves, with the reflection 

coefficient, Rq, varying from 0 to 1. a) fto vs Rq for 

0_=0°. b) fto vs Rq for 0„,u,=2O°. c) Ratio between 

estimated (fi) and true (E) incident wave energy vs Rq. d) 
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between the estimates and true values. Pluses and circles 
in panels c) cuid d) correspond to 0a,an=O° 20° 

simulations. 

Figure 5. Results of a model test for a reflection 
coefficient Ro=.5 and mean incidence angles varying from 

00-20°. The estimated incidence angle §nean is compared to 
^mran- The Solid line denotes perfect agreement. 

Figure 6. Estimated reflection coefficient :&o vs. frequency 
for all 18 data runs. The solid and dashed lines indicate 
low cuid high tide runs. 

Figure 7. Misfit e between and as a function of 
frequency for all 18 data runs. The solid euid dashed lines 
indicate low and high tide runs. 

Figure 8. The estimated phase lag $o between incident and 
reflected waves (relative to x=0) as a function of frequency 
for all 18 data runs. Circles and crosses indicate low and 
high tide runs. The solid line is a crude theoretical 
approximation of (|>o. 

Figure 9. The treuismission coefficient T as a function of 
frequency for all 18 data runs. The solid and dashed lines 
indicate low cuid high tide runs. 

Figure 10. The ratio of dissipated to incident energy flux 
as a function of frequency for all 18 data riins. The solid 
and dashed lines indicate low euid high tide runs. 

Figure 11. The ratios of reflected to incident (a), 
trauismitted to incident (b), and dissipated to incident (c) 
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Figur* 11. The ratios of reflected to incident (a), 
transmitted to incident (b), and dissipated to incident (c) 
energy flxixes vs. the incident energy flux for all 18 data 
rxins. The energy fluxes are integrated across the swell 
band (0.05-0.12). The circles and crosses represent low and 
high tide rvins. 
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Figure 2. 


27 





part OamO) real part 6nn(0) 





Flgur* 3. 
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